PROCEEDINGS 

OP SCIENCE 



Lessons learnt from INTEGRAL AGN 



Volker Beckmann", C. Ricci'', S. SoldK, J. Alfonso-Garzon'', T.J.-L. Courvoisier^ 
A. Domingo ', N. Gehrels^ P. Lubihski^, J.M. Mas-Hesse'', A.A. Zdziarski^ 

" Centre Frangois Arago, APC, IN2P3/CNRS, Universite Paris Diderot, 10 rue Alice Domon et 
Leonie Duquet, 75205 Paris Cedex 13, France 

^ ISDC Data Centre for Astrophysics, 16 Ch. d'Ecogia, 1290 Versoix, Switzerland 
'■' Laboratoire AIM - CNRS - CEA/DSM - Universite Paris Diderot (UMR 7158), CEA Saclay, 
DSM/IRFU/SAp, 91191 Gif-sur-Yvette, France 

Centro de Astrobiologia LAEX (CSIC-INTA), 28691 Yillanueva de la Canada, Madrid 
" Astrophysics Science Division, Code 661, NASA Goddard Space Flight Center, MD 20771, USA 
f Centrum Astronomiczne im. M. Kopernika, PL-00-716 Warszawa, Poland 



E-mail: beckmann@apc . univ-paris7 . f r 



The INTEGRAL mission provides a large data set for studying the hard X-ray properties of AGN 
and allows to test the unified scheme for AGN. We present results based on the analysis of 199 
AGN. A difference between the Seyfert types is detected in slightly flatter spectra with higher 
cut-off energies and lower luminosities for the more absorbed/type 2 AGN. When applying a 
Compton reflection model, the underlying continua (F ~ 1.95) appear the same in Seyfert 1 and 
2, and the reflection strength is /? ~ 1 in both cases, with differences in the inclination angle only. 
A difference is seen in the sense that Seyfert 1 are on average twice as luminous in hard X-rays 
than the Seyfert 2 galaxies. The unified model for Seyfert galaxies seems to hold, showing in 
hard X-rays that the central engine is the same in Seyfert 1 and 2 galaxies, seen under different 
inclination angle and absorption. Based on our knowledge of AGN from INTEGRAL data, we 
briefly outline open questions and investigations to answer them. In this context an ultra-deep 
12 Ms) extragalactic field can be a true legacy of the INTEGRAL mission in the area of AGN 
studies. 
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1. Introduction 

A longstanding discussion has been, whether Seyfert 1 and Seyfert 2 galaxies (or un-absorbed 
and absorbed Seyferts) indeed represent physically different types of objects, or if they can be 
unified under the assumption that they are intrinsically the same but seen from a different viewing 
angle with respect to absorbing material in the vicinity of the central engine, and that the difference 
in X-ray spectral slope can be explained solely by the absorption and reflection components. This 
unified model naturally explains the different Seyfert types in a way that the broad-line region is 
either visible (Seyfert 1) or hidden (Seyfert 2) possibly by the same material in the line of sight that 
is responsible for the absorption detectable at soft X-rays. On the other hand, the model has some 
problems explaining other aspects of AGN, for example, that some Seyfert galaxies change their 
type from 1 to 2 and back |[T],^. Also, the existence of Seyfert 2 galaxies that show no absorption 
in the soft X-rays, like NGC 3147 and NGC 4698 [|] cannot be explained by the unified model. 
The AGN surveys provided by SwiftlBKI [|] and INTEGRAL IBIS/ISGRI [|, |] have already led 
to the discovery that the fraction of absorbed and Compton thick sources is less than expected from 
cosmic X-ray background synthesis models (e.g. [Q]). With the ongoing INTEGRAL mission, it is 
now possible to compile a large sample of AGN for spectroscopic and correlation studies [§] and 
to probe the unified model for AGN. 

2. The intrinsic hard X-ray spectrum 

The effect that Seyfert 1 and low-absorbed objects appear to have steeper X-ray spectra than 
the Seyfert 2 and highly absorbed AGN was first noticed based on Ginga and CGRO/OSSE data 
[^] and later confirmed by various other studies. The difference in the hard X-ray spectral slope 
between Seyfert 1 and 2 has been a point of discussion ever since its discovery. 

A solution might be provided when considering the effects of Compton reflection on the hard 
X-ray spectrum. Recent analysis of a sample of 105 Seyfert galaxies using the spectra collected 
with BeppoSAX in the 2-200 keV band [ [To| ] provided no evidence of any spectral slope difference 
when applying more complex model fitting including a reflection component (PEXRAV). The mean 
photon index values found for Seyfert 1 and Seyfert 2 samples were F = 1.89 ±0.03 and F = 
1.80 ±0.05. The difference between types 1 and 2 is seen in this model in the different strength 
of the reflection component, with R = 1.2 ±0.1 and R = 0.9 ±0.1, and different cut-off energies 
of Ec = 230±22keV and Ec = 376±42keV, for Seyfert 1 and Seyfert 2, respectively. It has to 
be pointed out that spectral slope, reflection strength, and cut-off energy are closely linked. The 
IBIS/ISGRI data have a disadvantage over broad-band data when studying the spectral shape of 
the hard X-ray continuum, as they lack information about the spectrum below 1 8 keV. But on the 
positive side there is no influence by data taken at different epochs. Even when e.g. Swift/XRT 
and INTEGRAL data are taken simultaneous, the XRT data over only a very small fraction of the 
total IBIS/ISGRI exposure time, usually of the order of a few per cent. In those cases, where 
spectra are taken by more than one instrument at different times, the flux variability can mimic a 
stronger or weaker reflection component or cut-off energy. Although the spectra of AGN as seen 
by IBIS/ISGRI above 20 keV show little variability, significant changes in flux and spectral shape 
are indeed seen in the energy range below 10 keV. 
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Table 1: Results from spectral fitting of a simple power law, a cut-off power law, and a Compton reflection 
model (PEXRAV) to the stacked IBIS/ISGRI spectra of INTEGRAL AGN 
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The INTEGRAL data show consistent slopes for the spectra of unabsorbed / type 1 and ab- 
sorbed / type 2 objects already when a simple power-law model is used [^]. When applying the 
model used by Dadina (2007) to the stacked INTEGRAL spectrum of Seyfert galaxies, we get 
similar results: the underlying powerlaw appears to have consistent (within 2a) spectral slope for 
type 1 (r = 1.96) and type 2 (F = 1.91) objects and the same reflection strength 7? ~ 1.1, when 
applying different inclination angles of / ~ 30° and / ~ 60°, respectively. The data do not allow 
to determine the cut-off energy or inclination angle when fitting the reflection component. When 
fitting a simple cut-off power law, the INTEGRAL data show the same trend as the BeppoSAX sam- 
ple, i.e. a lower cut-off energy for Seyfert 1 {Ec = 86keV) than for Seyfert 2 {Ec = 184keV). 
It has to be taken into account, though, that the fit to the Seyfert 2 data is of bad quality, and 
that fixing the cut-off here to the same value as derived for the Seyfert 1, also leads to the same 
spectral slope. When fitting a reflection model to the stacked data, one gets a consistent pho- 
ton index of F ~ 1.95 and reflection strength /? ~ 1.3 for both absorbed and unabsorbed AGN. 
Table [l] summarises the results of the spectral fitting to the stacked ISGRI spectra of the vari- 
ous Seyfert subtypes. The values of F and R from the PEXRAV spectral modeling agree with 
the correlation R = (4.54 ± 1.15) x F - (7.41 ±4.51) found for the BeppoSAX AGN sample [0], 
which, for the INTEGRAL sample with F = 1.95, would lead to R = 1.4. This R{Y) correla- 
tion was first noted based on Ginga data for extragalactic and Galactic black holes, leading to 
/? = (1.4±1.2) X io-4r(i2-4±i.2) which in our case would result in a smaller expected reflec- 
tion component with R = 0.6 but within la of the value detected here. Absorbed and unabsorbed 
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sources show a consistent turnover at about Ec = lOOkeV when a cut-off power law model is 
applied. 

The observed dichotomy of different spectral slopes for type 1 and type 2 objects might there- 
fore be caused by data with too low significance, which do not allow to fit the reflection component, 
or in general by a strong dependence of the spectral slope on the choice of the fitted model. One 
aspect that has to be kept in mind is the dependence of the reflection strength R on the model ap- 
pUed and on the geometry assumed. Murphy & Yaqoob (2009) showed recently that their model of 
a reflection spectrum from a Compton-thick face-on torus that subtends a solid angle of 2n at the 
X-ray source is a factor of ~ 6 weaker than that expected from a Compton-thick disc as modelled 
in PEXRAV. Therefore, applying this torus model to the data presented here would result in a much 
larger relative reflection strength. 

A difference between type 1 and type 2 objects is seen in the average luminosity of these 
subclasses. For 60 absorbed Seyfert galaxies, the average luminosity is (L2o-i()okev) = 2.5 x 
10^^ergs^\ more than a factor of 2 lower than for the 74 unabsorbed Seyfert with redshift in- 
formation ((L20-iookev) = 6.3 X 10^^ crgs^'). The differences in luminosities are exactly the same 
when excluding the 5 Compton thick objects. 



3. Open questions 

INTEGRAL data have already led to a number of results concerning AGN at hardest X-rays. 
The local luminosity function seems to be well established now, the source population at this energy 
and spectral characteristics are now well described. Further observations, which might lead to a 
doubling of the number of INTEGRAL AGN will not change significantly the results in this field. 

There are several investigations that can be done exclusively by INTEGRAL. 

• as it has been shown in Tab. |lj, the cut-off energy and reflection strength cannot be fit si- 
multaneously when using the stacked spectrum of all INTEGRAL AGN. It also shows that 
the cut-offs in AGN spectra appear to be mainly at high energies, i.e. >80 keV. In order 
to constrain the cut-off, IBIS/ISGRI data on the brightest Seyferts are needed. This can be 
done using the brightest sources of the INTEGRAL AGN catalogue M] as it has been shown 



already in several investigations (e.g. [14, 15]). The INTEGRAL community should aim at 



completing pointed observations on the brightest AGN in order to form a complete sample. 

to measure the evolution of the INTEGRAL detectable AGN it is necessary to perform an 
ultra-deep (^ 12 Ms) extragalactic field observation. Several medium-deep surveys have 
already been performed, e.g. in the XMM-LSS field [|l^], in the Virgo region |17], and on 



the north ecliptic pole. INTEGRAL observation time should now be efficiently used in order 
to extend one of these into a ultra-deep field, which will be able to determine the evolution 
in redshift. Expectations for such a field range from 40 to 110 detectable AGN, depending 
on the assumptions and the evolutionary behaviour. This homogeneous sample would also 
be suitable to settle the debate about the true fraction of Compton thick sources among the 
hard X-ray detected AGN. 

detailed variability studies can be performed using a small well-defined sample of the bright- 
est AGN, which have sufficient long exposures and/or several observation periods and com- 
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billing the INTEGRAL data with those collected by other missions. Few studies like this have 



been performed, e.g. on 3C 273 and further work is in progress [ ]19| ] 



the question of spectral variability at hardest X-rays can also be addressed uniquely by IN- 
TEGRAL using the brightest detected AGNs. Such studies have been performed e.g. on 3C 
273 NGC 4151 [0]. Other bright objects, e.g. NGC 4945 and IC 4329A, are still to be 



investigated. In this context, short (200 ksec) follow-up observations on the ~ 10 brightest 
variable AGN would be useful in order to study their long term (on years time scale) spectral 
behaviour. 

INTEGRAL observations of blazars in outburst provide the highest quality hard X-ray spectra 



and help to constrain the spectral energy distribution [22, 23]. This is especially important 
as Fermi, AGILE, and TeV experiments rely heavily on simultaneous multiwavelength ob- 
servations. 

4. Conclusions 

The AGN population detected by INTEGRAL is dominated by Seyfert galaxies in the local 
((z) = 0.03) universe, with moderate X-ray luminosity ((L20-iookev) = 4 x 10^^ ergs"'). Taking 
into account the results on the Eddington ratios of the INTEGRAL detected AGN [||], Seyfert 1 
galaxies appear to have higher luminosities ((L2o-iookev) = lO^^ergs^') and Eddington ratio ((Asyi) 
0.064) than the Seyfert 2 galaxies ((L20-iookev) = 2.5 x lO'^^ergs-i, (Asy2) = 0.02). The underly- 
ing continuum of the hard X-ray spectrum appears to be consistent between different Seyfert types, 
both when a simple power-law model is applied and when considering the effects of Compton re- 
flection. Applying the PEXRAV reflection model with no high-energy cut-off, the Seyfert 1 and 2 
galaxies show the same underlying power law with F ~ 1 .95 and a reflection component of /? ~ 1 . 1. 
On the other hand, when applying a cut-off power law model to the stacked spectra, the Seyfert 1 
show lower cut-off energies {Ec = Se+j^keV) than the Seyfert 2 objects (Ec = 184+^2 keV). 
However, the bad quality of the fit in the latter case and the fact that fixing the cut-off to the 
value of the Seyfert 1 leads to a similar spectral slope might indicate that the spectra are in- 
trinsically indeed the same. The same differences as for different Seyfert classes are observable 
when considering the intrinsic absorption: the unabsorbed sources also have higher luminosi- 
ties ((L20-iookev) = 6.3 X lO'*^ crgs"') and Eddington ratio ((A) = 0.06) than the absorbed AGN 
((Z.20-iookev) = 2.5 X 10«ergs-i, (A) =0.015). 

The overall picture can be interpreted within the scenario of a unified model. The whole 
hard X-ray detected Seyfert population fills the parameter space of spectral shape, luminosity, and 
accretion rate smoothly, and only an overall tendency is seen in which more massive objects are 
more luminous, less absorbed, and accreting at higher Eddington ratio. More evidence for the 
unified scheme is that a fundamental plane can be found between the mass of the central object 
and optical and X-ray luminosity [^. The correlation takes the form Ly °^ L^-^M^^, similar to 



what is found in previous studies between radio luminosity Lr, Lx, and Mbh 1 24]. This links the 
accretion mechanism with the bulge of the host galaxy and with the mass of the central engine 
in the same way in all types of Seyfert galaxies. Further exploitation of the data is necessary to 
study the spectral characteristics of a bright sub-sample in a consistent way. In addition, variable 
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and bright Seyferts can be used to study spectral variability by peri'orming additional short (~ 
200 ks) observation with IBIS/ISGRI in the upcoming AOs. Evolutionary effects cannot be studied 
by INTEGRAL at the time being. Only if an ultra-deep 12 Ms) extragalactic field study is 
performed, we will be able to clarify how the hard X-ray AGN fraction evolves with redshift. Such 
a study would represent one of the important legacies of the INTEGRAL mission. 

References 

[I] Cohen, R. D., Puetter, R. C, Rudy, R. J., et al. 1986, ApJ 311, 135 

[2] Aretxaga, I., Joguet, B., Kunth, D., Melnick, J., & Terlevich, R. J. 1999, ApJ, 519, L123 

[3] Pappa, A., Georgantopoulos, I., Stewart, G. C, & Zezas, A. L. 2001, MNRAS, 326, 995 

[4] Tueller, J., Mushotzky, R. R, Barthelmy, S., et al. 2008, ApJ, 681, 113 

[5] Beckmann, V., Soldi, S., Shrader, C. R., Gehrels, N., & Produit, N. 2006, ApJ, 652, 126 

[6] Bassani, L., Landi, R., & MaUzia, A. 2007, ApJ, 669, LI 

[7] Treister, E. & Urry, C. M. 2005, ApJ, 630, 115 

[8] Beckmann, V., Soldi, S., Ricci, C, et al. 2009, A&A, 505, 417 

[9] Zdziarski, A. A., Johnson, W. N., Done, C, Smith, D., & McNaron-Brown, K. 1995, ApJ, 438, L63 
[10] Dadina, M. 2008, A&A, 485, 417 

[II] Dadina, M. 2007, A&A, 461, 1209 

[12] Zdziarski, A. A., Lubiriski, R, & Smith, D. A. 1999, MNRAS, 303, Lll 

[13] Murphy, K. D. & Yaqoob, T. 2009, MNRAS, 397, 1549 

[14] Beckmann, V., Courvoisier, T. J.-L., Gehrels, N., et al. 2008, A&A, 492, 93 

[15] Molina, M., Bassani, L., Malizia, A., et al. 2009, MNRAS, 399, 1293 

[16] Virani, S.N., Treister, E., Urry, CM., et al. 2010, in prep. 

[17] Paltani, S., Walter, R., McHardy, 1. M., et al. 2008, A&A, 485, 707 

[18] Soldi, S., Turler, M., Paltani, S., et al. 2008, A&A, 486, 41 1 

[19] Soldi, S., Ponti, G., Beckmann, V., & Lubinski, P. 2010, these proceedings, [arXiv : 1001 . 4348 ] 
[20] Chernyakova, M., Neronov, A., Courvoisier, T.J.-L., et al. 2007, A&A, 465, 147 
[21] Lubinski, P., Walter, R., Beckmann, V., et al. 2010, in prep. 
[22] Pian, E., Foschini, L., Beckmann, V., et al. 2006, A&A, 449, L21 

[23] Beckmann, V., Ricci, C, & Soldi, S. 2009, proc. of the Fermi Symposium 2009, 

[arXiv: 0912 .2254] 

[24] Merloni, A., Heinz, S. & Di Matteo, T. 2003, MNRAS, 345, 1057 



6 



